UT-722, 1995 



BRST Symmetric Formulation of a Theory 
^ ; with Gribov-type Copies 

+-> 

o 

Q \ Kazuo Fujikawa 

t> 

Department of Physics, University of Tokyo 
Bunkyo-ku, Tokyo 113, Japan 



in 

ON ■ Abstract 



23 



A path integral with BRST symmetry can be formulated by summing the 
Gribov-type copies in a very specific way if the functional correspondence between 



t and the gauge parameter u> defined by r(x) = f(A^) is "globally single val- 
•rH , 

^ ■ ued", where f(A^) = specifies the gauge condition. A soluble gauge model with 

h ; 

Gribov-type copies recently analyzed by Friedberg, Lee, Pang and Ren satisfies this 
criterion. A detailed BRST analysis of the soluble model proposed by the above au- 
thors is presented. The BRST symmetry, if it is consistently implemented, ensures 
the gauge independence of physical quantities. In particular, the vacuum (ground) 
state and the perturbative corrections to the ground state energy in the above model 
are analysed from a view point of BRST symmetry and i?^-gauge. Implications of 
the present analysis on some aspects of the Gribov problem in non-Abelian gauge 
theory, such as the 1/N expansion in QCD and also the dynamical instability of 
BRST symmetry, are briefly discussed. 
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1 Introduction 



The quantization of non-Abelian gauge theory is complicated by the presence of the 
so-called Gribov problem [4]. The Gribov problem in general suggests that the 
Coulomb gauge cannot completely fix the gauge due to the presence of more than one 
gauge configurations which satisfy the Coulomb gauge condition or , in certain circum- 
stances in compactified space-time, it even suggests the absence of the gauge configuration 
which satisfies the Coulomb gauge condition [2]. In Euclidean formulation of gauge the- 
ory, the Landau gauge also suffers from the same complications. Although the full details 
of the Gribov problem are not understood yet, a working prescription is figured out if 
one regards the Gribov problem as the appearance of several gauge copies( i.e., if one 
assumes that one can always find at least one gauge configuration which satisfies the 
Coulomb gauge condition). In fact, it was noted some time ago[§] that if the functional 
correspondence between r and the gauge parameter uj defined by 

r(x) = 0M£(x) (1.1) 

is "globally single valued", the path integral with BRST symmetry ||] is defined by sum- 
ming over Gribov copies in a very specific way. For a more general gauge fixing, (1.1) 
may be replaced by 

r(x) = f(A;(x)) (1.2) 

where f(A"(x)) = specifies the gauge condition. In (1.1) or (1.2) A^(x) stands for the 
gauge field obtained from A^x) by a gauge transformation specified by the gauge orbit 
parameter u(x). For an infinitesimal u>(x), one has 

A?{x) = Al{x) + d,co a (x) - gf abc Al{x)uj\x) (1.3) 

One may regard (1.1) or (1.2) as a functional correspondence between r and u 
parametrized by A^. The globally single valued correspondence between r and u, which 
is explained in more detail in Section 2, is required for any value of A^ to write a simple 
path integral formula. [ If one restricts oneself to only one of the Gribov copies by some 
means[0], one can also incorporate the idea of BRST symmetry. But a simple prescription 
which selectively picks up only one copy appears to be misssing at this moment.] 
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Recently a very detailed analysis of the Gribov problem was performed by Friedberg, 
Lee ,Pang and Ren|| on the basis of a soluble gauge model which exhibits Gribov-type 
copies. One of the main conclusions in || is that the singularity associated with the so- 
called Gribov horizon is immaterial at least in their soluble model and one may incorporate 
all the Gribov copies in a very specific way. They showed how this prescription works in 
the soluble model proposed by them. 

The present work is motivated by the fact that the model and the gauge choice in 
Ref.|| satisfy our criterion in (1.1) or (1.2). One can thus formulate a BRST invariant path 
integral for the model proposed in [8] by summing over Gribov-type copies ; the physical 
quantities thus calculated agree with those in Ref . || . Note that the BRST symmetry deals 
with an extended Hilbert space which contains indefinite metric in general, although the 
physical sector specified by BRST cohomology contains only positive metric. Since the 
BRST symmetry plays a fundamental role in modern gauge theory, we here present a 
detailed BRST analysis of the soluble model proposed in [||. 



2 BRST invariant path integral in the presence of 
Gribov copies 

In this Section we recapitulate the essence of the argument presented in ||. We start 
with the Faddeev- Popov formulation of the Feynman-type gauge condition P|[fLU[|. The 
vacuum-to-vacuum transition amplitude is defined by 

(+oo| - oo> = JvA;VC5(d»A; - C)A(A)exp{iS(A») - JL/ C(x) 2 dx} (2.1) 

where S(A^) stands for the action invariant under the Yang-Mills local gauge transforma- 
tion. The positive constant a is a gauge fixing parameter which specifies the Feynman-type 
gauge condition. [The equations in this Section are written in the Minkowski metric, but 
they should really be interpreted in the Euclidean metric to render the functional integral 
and the Gribov problem well-defined.] In the following we often suppress the internal 
symmetry indices, and instead we write the gauge parameter explicitly: A^ indicates the 
gauge field which is obtained from A^ by a gauge transformation specified by uj(x). The 
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determinant factor A (A) is defined by|| 

A(A)- 1 = JvuVC 5(d»A^-C)exp{-^Jc(x) 2 dx} 

« consti^ldeti-^-d^A^- 1 } (2.2) 

k ^ 

where the summation runs over all the gauge equivalent configurations satisfying d^A^f = 
, which were found by GribovfQ] and others [PJPIFJ]- Equation(2.2) is valid only for 
sufficiently small a, since the parameter A, C) defined by 

has a complicated branch structure for large C in the presence of Gribov ambiguities. 
Obviously the Feynman-type gauge formulation becomes even more involved than the 
Landau-type gauge condition. 

It was suggested in [|5| to replace equation (2.1) by 

(+oo| - oo) = iy VA^VC6(d^ - C)det[-^d^)exp{iS(A») -±-J C{xfdx} (2.4) 

The crucial difference between (2.1) and (2.4) is that (2.4) is local in the gauge space u(x) 
(i.e., the gauge fixing factor and the compensating factor are defined at the identical uo 
), whereas A (A) in (2.1) is gauge independent and involves a non-local factor in u as is 
shown in (2.2). As the determinant in (2.4) depends on A^, the entire integrand in (2.4) is 
in general no more degenerate with respect to gauge equivalent configurations even if the 
gauge fixing term itself may be degenerate for certain configurations. Another important 
point is that one takes the absolute values of determinant factors in (2.2) thanks to the 
definition of the 5-function, whereas just the determinant which can be negative as well 
as positive appears in (2.4). It is easy to see that (2.4) can be rewritten as 

(+oo| - oo) = i f VA^VBVcVc exp{iS{A^) + i f C g dx} (2.5) 



where 



C g = -d"B a A^ + id»c a (d^ - gf abc A b pc c + ^B a B a (2.6) 

with B a the Lagrangian multiplier field, and c a and c a the (hermitian) Faddeev-Popov 
ghost fields; f abc is the structure constant of the gauge group and g is the gauge coupling 
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constant. If one imposes the hermiticity of c a and c a , the phase factor of the determinant 
in (2.4) cannot be removed. The normalization constant N in (2.5) includes the effect of 
Gaussian integral over B in addition to N in (2.4), and in fact N is independent of a. 
Seeeq. (2.10). 

This C g as well as the starting gauge invariant Lagrangian are invariant under the 
BRST transformation defined by 



5 e c a = i6{g/2)f abc c b c c 
S e c a = 9B a 

5 e B a = (2.7) 

where 9 and the ghost variables c a (x) and c a (x) are elements of the Grassmann algebra , 
i.e., 9 2 = 0. This transformation can be confirmed to be nil-potent S 2 = 0, for example, 

M<W) = ° ( 2 - 8 ) 

One can also confirm that the path integral measure in (2.5) is invariant under (2.7). Note 
that the transformation (2.7) is "local" in the u parameter; precisely for this property, 
the prescription in (2.4) was chosen. 

To interprete the path integral measure in (2.4) as the path integral over all the gauge 
field configurations divided by the gauge volume, namely 

r VA^ 

(+oo| - oo> = / ? exp{iS(A")} (2.9) 

J gauge volume(u) M 

one needs to define the normalization factor in (2.4) by 

N = jvuVC5(d^A^-C)det[-^-d^A^}exp{-^-Jc(x) 2 dx} 
= Jvudet[^A-]exp{-^-J(d^A-) 2 dx} 

= Jvt exp{-^Jr(x) 2 dx} (2.10) 
where the function r is defined by 

T{x) = d»A"{x) (2.11) 



5 



and the determinant factor is regarded as a Jacobian for the change of variables from u(x) 
to t(x). Although we use the Feynman-type gauge fixing (2.11) as a typical example in 
this Section, one may replace (2.11) by 

r(x) = f(A%{x)) (2.12) 

to deal with a more general gauge condition 

f(A»(x)) = (2.13) 

It is crucial to establish that the normalization factor in (2.10) is independent of A^. 
Only in this case, (2.4) defines an acceptable vacuum transition amplitude. The Gribov 
ambiguity in the present case appears as a non-unique correspondence between t(x) and 
u>(x) in (2.11), as is schematically shown in Fig. 1 which includes 3 Gribov copies. The 
path integral in (2.10) is performed along the contour in Fig. 1. As the Gaussian function 
is regular at any finite point, the complicated contour in Fig. 1 gives rise to the same result 
in (2.10) as a contour corresponding to A^ = 0. In the present path integral formulation, 
the evaluation of the normalization factor in (2.10) is the only place where we explicitly 
encounter the multiple solutions of gauge fixing condition. [ If the normalization factor N 
should depend on gauge field A M , the factor N, which is gauge independent in the sense 
that we integrated over entire gauge orbit, needs to be taken inside the path integral in 
(2.4). In this case one looses the simplicity of the formula (2.4).] 

The basic assumption we have to make is therefore that (2.11) in the context of the 
path integral (2.10) is "globally single- valued" , in the sense that the asymptotic func- 
tional correspondence between u and r is little affected by a fixed A^ with d^A^ = 00; 
Fig.l satisfies this requirement. This assumption appears to be physically reasonable if 
the second derivative term of the gauge orbit parameter dominates the functional corre- 
spondence in (2.11), though it has not been established mathematically. To define the 
functional correspondence between u and r in (2.11), one needs in general some notion 
of norm such as L 2 -norm for which the Coulomb gauge vacuum is unique [3] [4]. The 
functional configurations which are square integrable however have zero measure in the 
path integral [1 1] , and this makes the precise analysis of (2.11) very complicated: At least 
what we need to do is to start with an expansion of a generic field variable in terms of 
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some complete orthonormal basis set (which means that the field is inside the L 2 -space) 
and then let each expansion coefficient vary from — oo to +oo(which means that the field 
is outside the L 2 -space). 

If the Gribov problem simply means the situation as is schematically shown in Fig. 
1, the prescription in (2.4) may be justified. The indefinite signature of the determinant 
factor in (2.4) is not a difficulty in the framework of indefinite metric field theory [ 12 1 [13 



since the determinant factor is associated with the Faddeev-Popov ghost fields and the 
BRST cohomology selects the positive definite physical space. On the other hand , the 
Gribov problem may aslo suggest that one cannot bring the relation(2.11) with fixed 
to d^A^ = by any gauge transformation 0. If this is the case, the asymptotic 
behavior of the mapping (2.11) is in general modified by A^ and our prescription cannot 
be justified. Consequently, the prescription (2.4) may be valid to the extent that one 
can always achieve the condition d^A^ = by means of suitable ( but not necessarily 
unique) gauge transformations. We also note that the limit a — > large in (2.5) and (2.10) 
corresponds to a very loose gauge fixing. 

3 BRST analysis of a soluble gauge model 

3.1, A SOLUBLE GAUGE MODEL 

The soluble gauge model of Friedberg, Lee, Pang and Ren|| is defined by 

C=±{[X(t)+gt{t)Y(t)]* + ^ (3.1) 

where X(t), for example, means the time derivative of X(t), and the potential U depends 
only on the combination X 2 + Y 2 . This Lagrangian is invariant under a local gauge 
transformation parametrized by u(t) , 

X u (t) = X(t) cos gu(t) -Y(t)smgu(t) 
Y u {t) = X(t)smgu(t) + Y(t) cos gu(t) 
Z"(t) = Z(t)+u(t) 

= m+"(t) (3-2) 
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The gauge condition ( an analogue of A = gauge ) 

C(t) = (3.3) 

or ( an analogue of A% = gauge ) 

Z(t) = (3.4) 

is well-defined without suffering from Gribov-type copies. However, it was shown in || 
that the gauge condition ( an analogue of the Coulomb gauge ) 

Z(t) - XX (t) = (3.5) 

with a constant A suffers from the Gribov- type complications. This is seen by using the 
notation in (3.2) as 

Z u (t)-XX u (t) = Z(t)+u(t)-\X(t)cosguj(t) + \Y(t)smguj(t) 

= u(t) + AVX 2 + F 2 [cos0(t) - cos(guj(t) + <f>(t))] = (3.6) 

where we used the relation(3.5) and 

X(t) = VX 2 + Y 2 cos cj)(t), Y(t) = VX 2 + Y 2 sin <j>{t) (3.7) 

From a view point of gauge fixing, u(t) = is a solution of (3.6) if (3.5) is satisfied. By 
analyzing the crossing points of two graphs in (u,rj) plane defined by 

f 

V = — iiu 
WX 2 + Y 2 

tj = cos(gu + (f>) — cos (f) (3.8) 

one can confirm that eq.(3.6) in general has more than one solutions for o>. 

From a view point of general gauge fixing procedure, we here regard the algebraic 
gauge fixing such as (3.3) and (3.4) well-defined; in the analysis of the Gribov problem 
in Ref. [2], the algebraic gauge fixing [3] is excluded. 

The authors in Ref.|| started with the Hamiltonian formulated in terms of the well- 
defined gauge = in (3.3) and then faithfully rewrote the Hamiltonian in terms of the 
variables defined by the "Coulomb gauge" in (3.5). By this way, the authors in |§ analyzed 
in detail the problem related to the Gribov copies and the so-called Gribov horizons where 
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the Faddeev-Popov determinant vanishes. They thus arrived at a prescription which sums 
over all the Gribov-type copies in a very specific way. As is clear from their derivation, 
their specification satisfies the unitarity and gauge independence. 

In the context of BRST invariant path integral discussed in Section 2, the crucial 
relation (2.11) becomes 

r(t) = Z"(t) - XX u (t) 

= w(t) +Z(t) - XX (t) cos gu(t) + XY(t) sin gu(t) (3.9) 

in the present model. For X = Y = , the functional correspondence between uj and r 
is one-to-one and monotonous for any fixed value of t . When one varies X(t), Y(t) and 
Z(t) continuously, one deforms this monotonous curve continuously. But the asymptotic 
correspondence between u(t) and r(t) at u(t) = ±oo for each value of t is still kept 
preserved, at least for any fixed X(t),Y(t) and Z(t) . This correspondence between oj(i) 
and r(i) thus satisfies our criterion discussed in connection with (2.11). The absence of 
terms which contain the derivatives of u(t) in (3.9) makes the functional correspondence 
in (3.9) well-defined and transparent. 

From a view point of gauge fixing in (3.6), this "globally single- valued" correspondence 
between u and r means that one always obtains an odd number of solutions for (3.6). The 
prescription in || is then viewed as a sum of all these solutions with signature factors 
specified by the signature of the Faddeev-Popov determinant 

det{-^-[Z"(t) - XX u (t)]} = det{[l + \gY"(t))8{t - 1')} (3.10) 



evaluated at the point of solutions , u = u^X 2 + Y 2 , (ft), of (3.6). The row and column 
indices of the matrix in (3.10) are specified by t and t', respectively. In the context of 
BRST invariant formulation, a pair-wise cancellation of Gribov-type copies takes place, 
except for one solution, in the calculation of the normalization factor in (2.10) or (3.21) 
below. 

3.2, BRST INVARIANT PATH INTEGRAL 

The relation (3.9) satisfies our criterion discussed in connection with (2.11). We can 
thus define an analogue of (2.5) for the Lagrangian (3.1) by 

(+oo| -oo) = i fdfi exp{iS{X u \Y" ', Z" \C) + * / ' C g dt} (3.11) 



where 

5(X",F",Z W ,D = J 'c{X u ,Y u ,Z?,?')dt 

= S(X,Y,Z,0 (3.12) 

in terms of the Lagrangian C in (3.1). The gauge fixing part of (3.11) is defined by 

C g = -PBC + B{Z" - AX W ) + pitc - ic(l + g\Y u )c + -B 2 (3.13) 

where a, (5 and A are numerical constants, and c and c are (hermitian) Faddeev-Popov 
ghost fields. B is a Lagrangian multiplier field. Note that C g is hermitian. The integral 
measure in (3.11) is given by 

dpi = VX u VY u VZ u V£ u VBVSDc (3.14) 

The Lagrangians C and C g and the path integral measure (3.14) are invariant under the 
BRST transformation defined by 



X^it^) 


= X u (t)-i6gc(t)Y u (t) 


Y u (t,e) 


= Y u {t)+i6gc(t)X u {t) 


Z w (t,B) 


= Z^^+iOcit) 




= ?(t)+i0c(t) 


c(t,0) 


= c(t) 


c(t,e) 


= c(t) + 9B(t) 



where the parameter 9 is a Grassmann number, 9 2 = 0. Note that 9 and ghost variables 
anti-commute. In (3.15) we used a BRST superfield notation: In this notation, the second 
component of a superfield proportional to 9 stands for the BRST transformed field of the 
first component. The second component is invariant under BRST transformation which 
ensures the nil-potency of the BRST charge. In the operator notation to be defined later, 
one can write , for example, 

X"(t, 9) = e-^X^it, 0)e eQ (3.16) 

with a nil-potent BRST charge Q, {Q, Q} + = 0. Namely, the BRST transformation is a 
translation in #-space, and 9Q is analogous to momentum operator. 
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In (3.1 1)~ (3.15), we explicitly wrote the gauge parameter ui to emphasize that BRST 
transformation is "local" in the u-space. For £ g in (3.13), the relation (3.9) is replaced 
by ( an analogue of the Landau gauge ) 



r(t) = + - XX u (t) 

= (3u{t) + u(t) + p£(t) + Z{t) - XX (t) cos gu(t) + XY(t) sin gu(t) (3.17) 

The functional correspondence between uj and r is monotonous and one-to-one for weak 
Z(t), X(t) and Y(t) fields; this is understood if one rewrites the relation (3.17) for 
Euclidean time t = —it E by neglecting weak fields as 



The Fourier transform of this relation gives a one-to-one monotonous correspondence 
between the Fourier coefficients of r and to for non-negative (3. The asymptotic functional 
correspondence between uj and r for weak field cases is preserved even for any fixed strong 
fields £(t), Z(t), X(t) and Y(t) for non-negative f3 to the extent that the term linear in 
uj(t) dominates the cosine and sine terms. The correspondence between r and u in (3.17) 
is quite complicated for finite u(t) due to the presence of the derivatives of ui{t). 

Thus (3.17) satisfies our criterion of BRST invariant path integral for any non-negative 
(3. The relation (3.9) is recovered if one sets /3 = in (3.17); the non-zero parameter 
(3 7^ however renders a canonical structure of the theory better-defined. For example, 
the kinetic term for ghost fields in (3.13) disappears for (3 = 0. In this respect the gauge 
(3.5) is also analogous to the unitary gauge. In the following we set (3 = a > in (3.13), 

C g = -aBC + B{Z W - AX") + aitc - ic(l + gXY^c + -B 2 (3.18) 

and let a — >• later. In the limit a = 0, one recovers the gauge condition (3.5) defined 
in Ref. ||. This procedure is analogous to i?^-gauge ( or the ^-limiting process of Lee 
and Yang|14| ) [p~5|j , where the (singular) unitary gauge is defined in the vanishing limit 



of the gauge parameter, £ — > 0: In (3.18) the parameter a plays the role of £ in R^- 
gauge. In i?g-gauge one can keep a ^ without spoiling gauge invariance, and the 
advantage of this approach with a ^ is that one can avoid the appearance of (operator 
ordering) correction terms || when one moves from the Hamiltonian formalism to the 
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Lagrangian formalism and vice versa. This point will be discussed in detail when we 
analyze perturbative corrections to the ground state energy in Section 4. 

By using the BRST invariance, one can show the A- independence of (3.11) as follows: 

(+oo|-oo) A+5A = (+oo|-oo) A -<5A-l fdti[B(t)X u (t)+igc(t)Y u (t)c(t)] exp{i f (C+C g )dt} 

TV J J 

(3.19) 

where we perturbatively expanded in the variation of C g for a change of the parameter 
X + 5X , 

£ g {\ + SX) = C g (X) - 5A[5(t)X w (t) + igc^Y^it)^)} (3.20) 
This expansion is justified since the normalization factor defined by (see eq.(2.10)) 

N = Jvt exp{-^-jT(x) 2 dt} (3.21) 

is independent of A provided that the global single- valuedness in (3.17) is satisfied. [See 
also the discussion related to eq.(5.1).] As was noted before, this path integral for N, 
which depends on a, is the only place where we explicitly encounter the Gribov-type 
copies in the present approach. By denoting the BRST transformed variables by prime, 
for example, 

X u( t y = X u,( t j _ ig Q c (t)Y"(t) (3.22) 
we have a BRST identity ( or Slavnov- Taylor identity [16]) 



1 

N 



dfic(t)X u (t) exp{i J(C + C g )dt} 
= 1 f dn'c(t)'X u (t)' exp{i ( (£' + C' g )dt} 
= iy dfJc(t)X u (t) exp{i J(C + C g )dt} 

+ fli f rf/i[ J B(t)X w (t) + igc(t)Y u (t)c(t)] exp{i f(C + C g )dt} (3.23) 

where the first equality holds since the path integral is independent of the naming of 
integration variables provided that the asymptotic behavior and the boundary conditions 
are not modified by the change of variables. The second equality in (3.23) holds due to 
the BRST invariance of the measure and the action 

dfi' = d/j,, 

C' + C' g = C + C g (3.24) 
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but 

c(t)'X u (t)' = c(t)X w (t) + #[5(t)X w (t) + igc(t)Y u (t)c(t)] (3.25) 
From (3.23) one concludes 

i f dn[B(t)X u (t) + ^c(t)F w (t)c(t)] exp{i /*(£ + C g )dt} = (3.26) 

and thus 

(+oo| - oo) a +5a = (+oo| - oo) A (3.27) 

in (3.19). This relation shows that the ground state energy is independent of the parameter 
A; in particular one can choose A = in evaluating the ground state energy, which leads 
to the gauge condition (3.4) without Gribov complications. 

In the path integral (3.11), one may impose periodic boundary conditions in time 
t on all the integration variables and let the time interval — > oo later so that BRST 
transformation (3.15) be consistent with the boundary conditions. 

The analysis in this sub-section is general but formal. In the next sub-section we 
convert the path integral (3.11) to an operator Hamiltonian formalism and analyze in 
detail the structure of the ground state and the gauge independence of physical energy 
spectrum. 

3.3, BRST ANALYSIS: OPERATOR HAMILTONIAN FORMALISM 

We start with the BRST invariant effective Lagrangian 



-e// 



= c + c g 

= \{[X«{t) + giityy^t)] 2 + \r>{t) - gC(t)X-(t)] 2 + [fr>(t) - C(t)] 2 } 



-U\(X"(t)f + (Y"(t)) 2 } 

-aB(t)C(t) + Bit^it) - \X"(t)) + ait(t)c(t) 

- ? c(t)[l + ^(t)]c(t) + |5(t) 2 (3.28) 

obtained from (3.1) and (3.18). A justification of (3.28), in particular its treatment of 
Gribov-type copies, rests on the path integral representation (3.11). In the following we 
suppress the suffix uj, which emphasizes that the BRST transformation is local in w-space. 
One can construct a Hamiltonian from (3.28) as 

H = hp* + p2 + Pl] + U(X 2 + Y 2 ) 
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+£G - B(Z - AX) + i-VcVc + ic(l + \gY)c - -aB 2 (3.29) 

a 2 

where 

Px = ^Ceff = X + gZY 

Py = ^ eff = Y-giX 
oY 

d 

Pz — — -jC P ff = Z — £ 

Pb = dB Ceff = ~^ 
P-c = -^C eff = iac 
d 

P c = -Qj£eff = -lac (3.30) 
and the Gauss operator is given by 

G = g(XP Y -YP x )+P z (3.31) 
We note that (p s y = —p s and (p c )t = —p c from (3.30). The quantization condition 

[Pb,B} = - (3.32) 



i 



implies 



% 



%B} = - (3.33) 



a 



and thus one may take the representation 

ia <9£ 



1 d 

B = — ^ (3.34) 



which is used later. 

We also note the quantization conditions 

% a 

{Pc,c} + = - -> {c,c}+ = - (3.35) 
% a 

The BRST charge is obtained from C e ff (3.28) via the Noether current as 

Q = C G- ip- c B (3.36) 
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The BRST charge Q is hermitian = Q and nil-potent 

{Q,Q} + = (3.37) 

by noting {c, c} + = {j9 g ,pc} + = {p s , c} + = 0. The BRST transformation (3.15) is gener- 
ated by Q, for example, 

e- eQ X(t)e eQ = X(t)-[6Q,X(t)} 
= X(t)-i6gc(t)Y(t), 
e- eQ c(t)e eQ = c(t)-[9Q,c(t)} 

= c(t) + 9B(t) (3.38) 

by noting 6 2 = 0. 

Some of the BRST invariant physical operators are given by 

X 2 + Y 2 , P 2 X + P*, 

L z = XP Y - YP X (3.39) 
The Hamiltonian in (3.29) is rewritten by using the BRST charge as 

Ol 

H = H + t{Q, & c } + + {Q, c(Z - XX} + + -{Q, Bc} + (3.40) 



with 



H = l -[P 2 x + Pi + Pf] + U(X 2 + Y 2 ) (3.41) 



One defines a physical state ^> as an element of BRST cohomology 

G fTer Q/Jm Q (3.42) 

namely 

Qty = (3.43) 

but \P is not written in a form \& = with a non-vanishing $. 
The time development of \I/ is dictated by Schroedinger equation 

i^!{t) = H*(t) (3.44) 
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and thus 

^ (St) = e- iH5t ^(0) 

= $(0)-iM$(0) 
= V(0)-i6tH V(0) 

ct 

-i5tQ{i£ Pc + c{Z - XX) + -Bc}^(0) 
~ *(0)-iM o $(0) (3.45) 

in the sense of BRST cohomology by noting (3.40) and Q^(0) = 0. Note that the 
Hamiltonian is BRST invariant 

[Q,H ] = [Q,H}=0 (3.46) 
If one solves the time independent Schroedinger equation 

H ^(0) = £tf(0) (3.47) 

with Q^(0) = 0, one obtains 



^(0)^-^(0) = 6-^*^(0)^(0) 

-iEt 



(3.48) 



by noting ^(0) f Q = 0. The eigen-value equation (3.47) is gauge independent and thus E 
is formally gauge independent, but the value of E is constrained by Q^(Q) = and thus 
we need a more detailed analysis. 

The basic task in the present BRST approach is to construct physical states \I/ satis- 
fying (3.42). We construct such physical states, in particular the ground state, as Fock 
states. For an explicit construction of physical states, we use the harmonic potential 
considered in Ref.131 



U(X 2 + Y 2 ) = y (X 2 + Y 2 ) (3.49) 



and we write H in (3.29) as 



H = \\Pl + P$] + ^r(X* + Yi) + l -( 9 L z Y 

H — p c p 5 + ice — \{Q, cX} + (3.50) 
a 
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with 



Pz = P z + gL z = G 
| = Z-gL z (3.51) 

Eq.(3.50) indicates that the freedom associated with Z has positive norm but the freedom 
£ has negative norm. We thus define 



X = 


1 / is 
- — (ax + a x ), 

yzuj 


Px = 


V 2 


Y = 


r—(a Y + a ] Y ), 


Py = 


V 2 


Z = 




Pz = 


Pz + gLz = ><J\[(b- 


1 = 




aB + Z = 




c = 


V2, (6 + 6l> ' 


Pc = 


-^-fc-c'l 


c = 




Pc = 





(3.52) 
with 

j //-J- (3.53) 
a V v 

In (3.52) the operator expansion of X, Y, c, and c is the standard one, but the expansion 
of Pz i and aB + Z is somewhat unconventional. 

The canonical commutators are satisfied by postulating 

[a x ,a x ] = 1, 

[ay, a},] = 1, 
17 



[d,rf] = 1, 

[b,b ] ] = -1, 
{c,c}+ = -i, 

{ct,g} + = i (3.54) 



and all other commutators are vanishing. Here we defined 



/2u 

so that 



b = b-i^=gL z (3.55) 



^ = (-0^(5-^), 

afi + Z = -— (6 + ftt) 
V2z/ 



The operator 6 also satisfies 

[6, ftt] = _i (3.56) 

The variables b in (3.54) and b in (3.56) carry negative norm. The Hamiltonian (3.50) is 
then written as 

H = u(a x a x + a\a Y + 1) + ^{gL z ) 2 

+v[d ] d - tfb + 1] + iv[c ] c - c j c + i] 

+\{Q,cX} + (3.57) 

with the BRST charge 

Q = cG — ip 5 B 

= iv[c ] (b - d) - c(b - d)t] 

= iv[$(b-d)-c(b-d)i] + yl^(& + c)gL z (3.58) 

and 

L z = XP Y - YP X = i(a x a Y - a Y a x ) (3.59) 
The BRST charge in (3.58) is nil-potent 

Q 2 = u 2 c ] c[b -d,tf- d f ] = (3.60) 
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by noting (3.54) and (3.56). 

We thus define the (physical) ground state at t = by 



6|0) = d|0)=0, 
c|0) = c|0) = 0, 
a x \0) = a Y \0) = (3.61) 

which ensures 

Q|0)=0 (3.62) 

The zero-point energy of Z and £ and the zero-point energy of c and c in (3.57) cancel 
each other for the state |0). When one defines a unitary transformation [8] 

a x = -^(a x -ia Y ) 

a Y = —j={—ia x + a Y ) (3.63) 
v2 

we can write the physical part of H in (3.57) as 

H phys = u(a x a x + a Y a Y + 1) + ^{gL z ) 2 

= ui{d x d x + a Y a Y + 1) + — (a x a x — a Y a Y ) 2 (3.64) 
If one recalls the relation 

1 d 

V2z/ lot, 

b-d = -l=[-i(P z + gL z ) + v^-] (3.65) 
V Iv loq, 

where we used aB = -Jj^ in (3.34), the state |0) in (3.61) is required to satisfy 



{^-e + i^}|0>z € = 0, 

{ ^ + 7^ + ^ }|0) ^ = ° (3 ' 66) 



Eq.(3.66) has a solution 



= N e W {±-{i-iuZ) 2 -igL z Z} (3.67) 
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which depends on Lz ; Nq is a normalization constant. The inner product of this state 
needs to be defined by means of a 90 degree rotation in the variable £ 

(0, L z \0, L z ) = iV 2 J dZd^ e ' uZ2 (3.68) 

which reflects the fact that the ^-variable carries negative norm. If one projects the state 
|0, Lz)zs, to the one with p^ = by Fourier transformation, which is the general procedure 
of Dirac in his treatment of singular Lagrangian [17], one obtains 

J dt%Lz)zi~er iaLzZ (3.69) 

This is the physical ground state naively expected for (Z, P z ) system for a given L z on 
the basis of invariance under the Gauss operator in (3.31), and it is independent of the 
gauge parameter v = y^. 

Other sectors of the ground state \0)xyec are constructed in a standard manner 

d X \0)xYcc = d Y \0)xYcc = c\0)xYcc = c|0) XYcc = (3.70) 

and the entire ground state in (3.62) is written as 

|0> = \0,L z )zt® \0)xycc (3.71) 

where Lz = a x ax — SySy is replaced by a c-number by acting Lz on \0)xycc- In the 
present case, L z = 0. 

The ground state thus defined has a time dependence described by H in (3.57) 

e - iH5t \0) = |0) -iu5t\0) + 8t\QcX\0) 

~ \0)-iw5t\0) (3.72) 

in the sense of BRST cohomology by noting (3.62). Thus we can represent the vacuum 
(or ground) state in Ker Q/Im Q by |0) in (3.71) for any time, and we have 

(0\ e - lHt \0) = e- lut (3.73) 

Excited physical states are represented by 

' '4) ni («y) n2 10) = |0, L z = rai - n 2 ) Z( ® ^i=(4) ni (4)" 2 |0)xy 5c (3.74) 
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One obtains the eigenvalue of H in (3.57) and (3.64) as 



+\QcX^J==(tf x r(tf Y r\0) 
\Jn\\n 2 \ 

~ Mm + n 2 + l) + £(m - n 2 ) 2 ]^i=(4) ni (4) n2 |0) 

(3.75) 

in the sense of BRST cohomology, since 

Q^i=(4r(4r|0) = (3.76) 
This (3.76) is confirmed by using the expression of Q in (3.58) as 
„-e Q 1 ^) ni (4) n2 |0> 



1 

\ /n V n 2 ] - 

1 '^) ni (4-r|0> (3.77) 



= e 



where we used 



L z = a x a x - a Y a Y , 



b\0,L z = ni - n 2 )zi = {b-i—7==g(n 1 -n 2 )}\0,Lz^n 1 -n 2 )z^ 



= (3.78) 

The second relation in (3.78) is regarded as a constraint on the ground state of b for a given 
value of L z , which is a manifestation of the Gauss constraint in the present formulation. 

The states which include b\ d\ c^and c excitations become unphysical and are re- 
moved by BRST cohomology. For example, the state 

-v{b-dy\0) = Q"c j \0) (3.79) 

is BRST invariant and has energy (for A = ) 

HQ$\0) = (lo + u)Qc j \0} (3.80) 
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but it is obviously excluded by BRST cohomology. The unphysical excitations 
b\ d\c\ and form the components of non-trivial BRST superfields Z u (t,9) and c(t,9) 
in (3.15), 

Z^{t,9) = Z"(t)+i9c(t) 

- — {d + d ] )+i9^(c + c ] ) 



c(t,9) = c(t) + 6B(t) 

= ^(Z+$) + 9v^[b-d + (b-dy} (3.81) 

A characteristic property of these non-trivial superfields is that the second components of 
the superfields, which are BRST transform of the first components, contain terms linear 
in the elementary field. The basic theorem of BRST symmetry is that any BRST invariant 
state which contains those unphysical degrees of freedom ,b\d),& and c , is written in 
the BRST exact form such as in (3.79) and thus it is removed by BRST cohomology. See 
Ref.[13]. 

The present BRST analysis is in accord with an explicit construction of physical states 
in Ref.p]. One can safely take the limit a — > (or v = ^= — > oo) in the physical sector, 
though unphysical excitations such as in (3.80) acquire infinite excitation energy in this 



limit just like unphysical excitations in gauge theory defined by R? -gauge |15 



4 Perturbative calculation in path integral 

It has been shown in M that the correction terms arising from operator ordering plays 
a crucial role in the evaluation of perturbative corrections to ground state energy in 
Lagrangian path integral formula. This problem is often treated casually in conventional 
perturbative calculations; a general belief ( and hope) is that Lorentz invariance and 
BRST invariance somehow takes care of the operator ordering problem. In the following, 
we show that BRST invariance and T*-product prescription reproduce the correct result 
of Ref . || provided that one uses a canonically well-defined gauge such as .R^-gauge with 
a 7^ in (3.18). This check is important to establish the equivalence of (3.11) to the path 
integral formula in Ref.||. See also Ref. [18]. If one starts with a = from the on-set, 
one needs correction terms calculated in Ref. ||. 
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To be precise, what we want to evaluate is eq.(3.11), namely 




(4.1) 



We define the path integral for a sufficiently large time interval 



t E [T/2, -T/2] 



(4.2) 



and let T — > oo later. In the actual calculation, there appear two important aspects which 
need to be taken into account: 

(i) We impose periodic boundary conditions on all the variables so that BRST transfor- 
mation (3.15) is well-defined including the boundary conditions. 

(ii) In the actual evaluation of the path integral as well as Feynman diagrams , we may 
apply the Wick rotation and perform Euclidean calculations. 

The exact ground state energy of (4.1) is given by eq.(3.73) as 



Namely, we have no correction depending on the gauge parameter A and the coupling 
constant g. As was already shown in (3.27), the absence of A dependence is a result of 
BRST symmetry. This property is thus more general and, in fact, it holds for all the 
energy spectrum of physical states; this can be shown by using the Schwinger's action 
principle [18] and the definition of physical states in (3.42). The perturbative check of A- 
independence or Slavnov- Taylor identities in general is carried out in the standard manner. 
On the other hand, the absence of g- dependence is an effect of more dynamical origin. 
In this Section we concentrate on the evaluation of (^-dependence by taking a view that 
the A- independence has been generally established in (3.27). 

We first perform the Gaussian path integral over Invariable by noting 



E = u 



(4.3) 



C. 



■gauge 



-B 2 + aB£ + B(Z -AX) 

fu 1 1 

-[B + -« + Z- AX)] 2 -—[at + z- AX] 2 
la la 

-^K + Z-AX] 2 



(4.4) 
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We thus consider the path integral 

(+00 1 - 00) = VXVYVZViVcDc exp{i J £ eff dt} (4.5) 
where N is the original normalization constant in (3.21), and 

£*ff = \{m)+gmyit)] 2 + [Y{t)- g mxm-^[x{tf+Y{t)^ 
+\z{tf - ±[z(t) - Axw] 2 - |ew 2 + \m 2 

+\i(t)X(t) + ait(t)c(t) - ic{t){\ + g\Y(t))c(t) 
= C + Cj (4.6) 

with 

C„ = \[X(tf + Y(t) 2 ]-^{Xhf + Y(t) 2 ] 

+aic(t)c(t) - ic(t)c(t), (4.7) 
Cj = gmiXW(t)-Y(t)X(t)] + ^atflX(tf + Y(tf] (4.8) 

where we set A = in the final expressions in (4.7) and (4.8). 

The propagators are defined by Co in (4.7) in the standard manner as 

(T*X(ti)X(t 2 )) = (T^Y^Y^)) 

dk ,- t 



2tt 

2vr 
' dk 



fa) 




k 2 


— cj 2 + ie 




i 


k 2 


— 1/a + ie 


t 2 ) 


—i/a 


A; 2 


— 1/a + ie 


fa) 


-1/a 



where we took T — > 00 limit in the evaluation of those propagators; this procedure is jus- 
tified for the evaluation of corrections to the ground state energy since all the momentum 
integrations in (4.11) below are well-convergent. 

Up to the second order of perturbation in £7, we have 

(+oo| - 00) = e */^*{l + i J dt{T* l -g^(tf[X(tf + Y{tf]) 
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+ 



(l) 2 



J dhdhiT+g^ihKih^XitjYitJ-YitjXih 



2! 

x[X(t 2 )Y(t 2 ) -Y(t 2 )X{t 2 )})} 
= jffdli e*J c ° dt {l + tg 2 J dt(T*Z(t)Z(t))(T*X(t)X(t)) 

+ (i)V/ *i^ 2 <^(ti)e(*2))<r*x(* 1 )x(t 2 ))(T*y(* 1 )y(t 2 )) 

- «V / *icft2<r*e(*i)e(*2)><r*x(t 1 )x(t 2 )>(r*y(t 1 )y(t 2 )>K4.io) 



by taking the symmetry in X and y into account. 

The terms of order g 2 in (4.10) are written by using the propagators in (4.9) as 

dk 1 r dl 1 



. i 2 . r r dk 1 r 
{--g )J dtj 7^ k 2 + 1 / a J 2nP + uj 2 

, % o, r , r dk 1 /■ I 2 

+i a g) J dt J^WTT^J 



2tt k 2 + 1/a J 2nl 2 + uj 2 (k + I) 2 + lu 2 
dk 1 r dl 1 1(1 + fc) 
2^k 2 + l/aJ 2^l 2 + u 2 (k + l) 2 + u 2 
dk 1 r dl 1 



27rA; 2 + l/a7 2tt Z 2 + cj 2 
d/c 1 r dl 1 



2tt A; 2 + 1/a 7 2vr Z 2 + cj 2 
i 2 r dk 1 r dl u 2 1 
+ J 2^k 2 + l/aJ 2^l 2 + u 2 {k + l) 2 + u 2 

+ ( V)T / - 1 / - — * ( * + fc) (4.11) 

y a } J 2irk 2 + l/aJ 2tt I 2 + u 2 {k + I) 2 + u 2 1 ; 

after the Wick rotation. We here note that the T*-product prescription is crucial in 

obtaining (4.11) ; the T*-product commutes with the time derivative operation, which 

is intuitively understood from the fact that the basic path integration variables are field 

variables ( and not their time derivatives) in the Lagrangian path integral formula [19]. 

In this approach, the conventional T-product is defined from T*-product via the Bjorken- 

Johnson-Low prescription [20] . 

The first two terms in (4.11) cancel each other. The last two terms in (4.11) can be 

evaluated as follows: In the third term in (4.11), one can rewrite the integrand as 

UP 2 1 1 2 -1 1 -1 1 

I 2 + uo 2 (k + I) 2 + uo 2 ~ { 2i } { l + iiu + l-iu )[ l + k + itu + 1 + k-iJ 

2% I + iuj I + k — iu I — iul + k + ico 

= (-) 2 \— - + - - ) (4 12) 

K 2 J l l + icol + k-iuj l-iul + k + iuj i v 1 
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since the poles displaced in the same side of the real axis do not contribute to the integral 
over / dl; the contour can be shrunk to zero without encircling poles. Similarly, in the 
last term in (4.11), 

1 1(1 + k) 1 1 1 1 1 



Z 2 + cj 2 (fc + I) 2 + cj 2 v 2 7 v / + icu l-iuo' y l + k + iuj l + k-iu' 

^ (-) 2 \ — + — } (4.13) 

v 2 7 l l + iuul + k-iuu I - iool + k + iu v ' 

The last two terms in (4.11) thus cancel each other, and we obtain only the lowest order 
contribution in (4.11) 

(+oo|-oo> = ^Jduc^ c ° dt 

1 1 det[ad 2 + 1] 



= const x 



[2i sin(^T/2)] 2 N ^det[8 2 t + l/a]det[ad 2 t + 1] 



= const x — — ; — for T — > oo, (4.14) 

[2isin(cuT/2)] 2 J v ; 

where the determinant factors coming from Z, ^, c and c integration combined with the 
normalization constant iV cancel completely among themselves. The last expression in 
(4.14) is a standard path integral of harmonic oscillators X and Y with periodic boundary 
conditions [21], and it may be expanded as 

1 oo 

const x — - = const x e~ wT (l + V e"^ T{ni+n2) ) (4.15) 

[2.sin(^T/2)] 2 1 n ± =0 ] 1 J 

for T — > oo, where the summation over the no n- negative integers rii and n 2 excludes the 
case ni = ri2 = 0. To be precise, the T-dependence of the normalization constant iV needs 
to be taken into account to obtain the last expression of (4.14) [21]. The ground state 
energy is then obtained from 

(+oo|-oo) = lim(0|e- iH [ T / 2 -(- r / 2 )]|0> 

= lim const x e~ lujT (4.16) 

T— >oo 

which is justified for T = — iT E and T E — > oo in Euclidean theory. We thus obtain the 
ground state energy 

E = u (4.17) 

to be consistent with (3.73). 
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5 Discussion and conclusion 



The BRST symmetry plays a central role in modern gauge theory and the BRST invariant 
path integral can be formulated by summing over all the Gribov-type copies in a very 
specific manner provided that the crucial correspondence in (2.11) or (3.17) is globally 
single valued [5]. This criterion is satisfied by the soluble gauge model proposed in Ref. [8], 
and it is encouraging that the BRST invariant prescription is in accord with the canonical 
analysis of the soluble gauge model in Ref. 0. The detailed explicit analysis in Ref.|| and 
the present somewhat formal BRST analysis are complementary to each other. In Ref. M, 
the problem related to the so-called Gribov horizon , in particular the possible singularity 
associated with it, has been analyzed in greater detail; this is crucial for the analysis of 
more general situation. On the other hand, an advantage of the BRST analysis is that one 
can clearly see the gauge independence of physical quantities such as the energy spectrum 
result of BRST identity. 

The BRST approach allows a transparent treatment of general class of gauge con- 
ditions implemented by (3.18). This gauge condition with a ^ renders the canonical 
structure better-defined, and it allows simpler perturbative treatments of the problems 
such as the corrections to the ground state energy. Our calculation vis-a-vis the explicit 
canonical analysis in Ref. |J may provide a (partial) justification of conventional covariant 
perturbation theory in gauge theory, which is based on Lorentz invariance( or T*-product 
) and BRST invariance without the operator ordering terms. 

Motivated by the observation in Ref. [8] to the effect that the Gribov horizons are not 
really singular in quantum mechanical sense, which is in accord with our path integral 
in (3.21), we would like to make a speculative comment on the role of Gribov copies in 
QCD. First of all, the topological phenomena related to instantons for which the Coulomb 
gauge is generally singular may be analyzed in the temporal gauge A = 0; this gauge is 
relatively free from the Gribov complications [3]. A semi-classical treatment of instantons 
with small quantum fluctuations around them will presumably give qualitatively reliable 
estimates of topological effects. 

As for other non-perturbative effects such as quark confinement and hadronic spec- 
trum, one may follow the argument of Witten [22] on the basis of 1/N expansion in 
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QCD[23]; he argues that the 1/N expansion scheme comes closer to the real QCD than 
the instanton analysis in the study of hadron spectrum. If this is the case, one can analyze 
the qualitative aspects of hadron spectrum on the basis of a sum of (an infinite number of) 
Feynman diagrams. This diagramatic approach or an analytical treatment equivalent to it 
in the Feynman-type gauge deals with topologically trivial gauge fields but may still suffer 
from the Gribov copies, as is suggested by the analysis in Ref. [4]: If one assumes that 
the vacuum is unique in this case as is the case in L 2 -space, the global single- valuedness 
in (2.11) in the context of path integral (2.10) will be preserved for infinitesimally small 
fields Ap. By a continuity argument, a smooth deformation of A^ in L 2 -space ( or its 
extension as explained in Section 2) will presumably keep the integral (2.10) unchanged. 
If this argument should be valid, our path integral formula in (2.4) would be justified. If 
our speculation is correct, the formal path integral formula (2.5) will provide a basis for 
the analysis of some non-perturbative aspects of QCD. 

On the other hand, the Gribov problem may also suggest the presence of some field 
configurations which do not satisfy any given gauge condition in four dimensional non- 
Abelian gauge theory [2]. For example, one may not be able to find any gauge parameter 
uj(x) which satisfies 

&>A»(x) = (5.1) 

for some fields A^. Although the measure of such field configurations in path integral is 
not known, the presence of such filed configurations would certainly modify the asymp- 
totic correspondence in (2.11). In the context of BRST symmetry, the Gribov problem 
may then induce complicated phenomena such as the dynamical instability of BRST sym- 
metry [24]. If the dynamical instability of BRST symmetry should take place, the relation 
corresponding to (3.27), which is a result of the BRST invariance of the vacuum, would 
no longer be derived. In the framework of path integral, this failure of (3.27) would be 
recognized as the failure of the expansion (3.19) since the normalization factor N in (3.21) 
would generally depend on not only field variables but also A if the global single- valuedness 
in (3.17) should be violated. 

Finally, we note that the lattice gauge theory [25], which is based on compactified field 
variables, is expected to change the scope and character of the Gribov problem completely. 
The Gribov problem is intricately related to the difficult issue of the non-perturbative 
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continuum limit of lattice gauge theory. [ Note that the perturbative continuum limit is 
not quite relevant in the present context, since we do not worry about the Gribov problem 
much in weak field perturbation theory]. 
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